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Abstract 

We propose that dark energy in the form of a scalar field could effectively 
couple to dark matter inhomogeneities. Through this coupling energy could be 
transferred to/from the scalar field, which could possibly enter an accelerated 
regime. Though phenomenological, this scenario is interesting as it provides 
a natural trigger for the onset of the acceleration of the universe, since dark 
energy starts driving the expansion of the universe when matter inhomogeneities 
become sufficiently strong. Here we study a possible realization of this idea by 
coupling dark energy to dark matter via the linear growth function of matter 
perturbations. The numerical results show that it is indeed possible to obtain a 
viable cosmology with the expected series of radiation, matter and dark-energy 
dominated eras. In particular, the current density of dark energy is given by the 
value of the coupling parameters rather than by very special initial conditions for 
the scalar field. In other words, this model - unlike standard models of cosmic 
late acceleration - does not suffer from the so-called “coincidence problem” and 
its related fine tuning of initial conditions. 

Keywords: Cosmology; Dark Energy; Dark Matter; Large-Scale Structure of 
the Universe. 


1. Introduction 

According to the standard model of cosmology - confirmed by the latest 
results from ESA’s Planck mission [1] - about 5% of the energy content of 
the universe is made of ordinary baryons, that is, of particles belonging to the 
standard model of particle physics, recently glorified by the discovery of the 
Higgs boson [2, 3]. The dark sector accounts for the remaining 95%. More 
precisely, roughly 25% consists of a yet-undetected matter component, which is 
thought to be a massive particle of non-baryonic nature that interacts through 
weak interaction and gravity only. It is dubbed “cold dark matter”. Finally, 
dark energy is responsible for the missing 70% of the energy content. The best 
candidate for dark energy to date is the so-called “cosmological constant”, which 
is basically the energy of the vacuum and, in general relativity, is an arbitrary 
constant of nature. Its fundamental property - gravitational repulsion - causes 
the acceleration of the expansion of the universe [4, 5]. 

Faced with the formidable challenge of accounting for not only one but two 
unknown components, cosmologists studied at great depth dynamical models 
of dark energy [6, 7] and a possible interaction between dark energy and the 
other fields [see 8, 9, 10, and references therein], hoping to shed light on the 
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nature of the dark sector. In particular, a coupling between dark energy and 
dark matter [11, 12, 13, 14] is suggested by the fact that the individual energy 
densities of the two components are today of the same order of magnitude. 
The latter is the so-called “coincidence problem”, as an incredible fine tuning 
of the value of the cosmological constant is necessary for the standard model 
to explain observations. However, a careful dynamical analysis of coupled dark 
energy models has shown that a fine tuning of the dark energy potential is 
unavoidable in order to have a viable cosmology [see, for example, 9]. 

Here we propose that dark energy could effectively couple to dark matter 
inhomogeneities. Indeed, if this is the case, the dark energy evolution will 
be altered - possibly causing the universe acceleration - at late time when 
inhomogeneities become strong. In other words, the trigger for acceleration 
is not given by tuned initial conditions but rather by the natural evolution of 
inhomogeneities. We then take a phenomenological approach and consider a 
possible realization of this general idea by coupling dark energy to dark matter 
via the linear growth function of matter perturbations. Our aim is to understand 
if such interaction can indeed give a viable dark energy dynamics free of fine 
tuning. 

The outline of this paper is as follows. First we introduce the model in 
sections 2-4. Next we discuss the limiting cases in section 5 and the evolution 
in section 6 . Conclusions are given in section 7. 

2. Model 

Within General Relativity the total energy-momentum tensor is conserved 
and a possible interaction between dark matter and a scalar field can be modeled 
through an interaction current which transfers energy and momentum from 
one source to the other: 

V a T^ = Q? V a T f = -QP . (1) 

In the individual conservation equations the interaction current has opposite 
sign so that the total energy-momentum tensor is conserved. 

From equation (1) (Q® = 0 because of homogeneity and isotropy) one can de¬ 
rive the following system of dynamical field equations describing the background 
evolution of the well known Coupled Dark Energy model [see, e.g., 9]: 


2 8nG 

H — ^ \Pm + P<t> + Pr) 5 

(2) 

H + H — (p m + + 3 + 2 p r ), 

(3) 

p<l> + 3 H{p</> + pj,) = — Q° , 

(4) 

Pm + 3 Hp m = Q° , 

(5) 

p r -h 4Hp r = 0 . 

(6) 


Curvature and the (subdominant) baryons are neglected, an overdot represents 
a derivative with respect to the cosmic time t, H = a/a is the Hubble function, 
and Mpi = l/y/8nG is the reduced Planck mass with G the Newton constant 
(we set c = h = 1). The first two equations are, respectively, the Friedmann 
and the acceleration equation, while the last three equations are the continuity 
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equations for the scalar field, matter and radiation, respectively. The scalar 
field energy density, pressure and equation of state are, respectively: 

P<t>=\<t> 2 + v (4'), P4>=\(f 1 = —■ ( 7 ) 

Z Z P(j) 


3. Interaction 

In the literature two classes of interactions - with several variations - have 
been studied: 

Q? = QT m V^/Mpi, (8) 

Qu = Q T m H , (9) 

where T m = —p m is the trace of the dark-matter energy-momentum tensor so 
that in equations (4) and (5) one has: 

Q°i = Q Pm HMpi , ( 10 ) 

Qu = Q Pm H . (11) 

The first class follows from the interaction Lagrangian = m((t>)ip'ip with 
m{(j >) = moe^/ Afp1 , 1 that is, the mass of the dark matter field ip is a function of 
the scalar field (p, [see 11, for details]. The second class is instead an arbitrary 
parametrization, and equation (1) is seen as a phenomenological description for 
an effective interaction [14]. As said earlier, careful dynamical analyses have 
shown that fine tuning of initial conditions is unavoidable in order to have a 
viable cosmology [see, for example, 9]. 

Motivated by the fact that the evolution of inhomogeneities could trigger at 
the right time the onset of a dark-energy dominated universe, we now discuss 
the possibility that the dark-energy scalar field effectively couples to matter in¬ 
homogeneities. 2 The coupling is realized by considering an interaction like (10), 
in which the coupling is a function of the linear growth function of density 
perturbations 5 m , Q = so that: 

Q° = is6"p m <p/M P1 . (12) 

The interaction (12) couples equations (4-5) to the linear perturbation equa¬ 
tion (13). The use of S m is perhaps the simplest way to quantify the overall 
growth of inhomogeneities in the universe. The dimensionless coupling param¬ 
eter v sets the strength of this interaction. The exponent n parametrizes the 
dependence of the effective coupling parameter Q = v 6^ on the linear growth 
function S m . For the sub-horizon scales at which a canonical scalar field remains 
homogeneous, the growth function does not depend on the wavenumber k and 
can be obtained numerically from the standard equation [9]: 

5 m + 2 HS m - AirGpmSm = 0 . (13) 


1 Generally, the coupling is Q = d In 'rn(p) / do. 

2 The idea that the evolution of inhomogeneities could trigger the acceleration of the uni¬ 
verse has been discussed in the literature regarding the so-called backreaction proposal, ac¬ 
cording to which late-time inhomogeneities could affect the average expansion rate of the 
universe, possibly explaining away dark energy (see, e.g., the special focus issue [15]). 
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As usual, we normalize 8 m to unity at the present time, S m (to) = 1. It is inter¬ 
esting to note that the coupling parameter v does not appear as the correction 
due to v is of order n +1 in 8 m . Therefore, the growth of perturbation is altered 
only by the different background evolution, meaning that this model could agree 
with growth data for arbitrarily large values of v. 

3.1. Discussion 

Although in the following we take the phenomenological approach in which 
we assume the validity of (12) and examine its cosmological implications, it is 
interesting to discuss which mechanism could be behind such an interaction. 

Physical processes that are linked to the growth of S m are necessarily nonlin¬ 
ear, such as collapse and shell crossing. However, quantities associated to scalar 
perturbations do not seem suitable for the coupling (12). Indeed, local values 
of e.g. matter density and velocity divergence can be arbitrary high in the past 
even if S m was initially very small. The obvious quantity that at early times 
is absent (or negligible) is the vorticity (or curl) of the velocity field, 3 which is 
generated at very non-linear stages of the dark matter collapse. More precisely, 
it is produced when, at shell crossing, the single-stream irrotational description 
of the dark matter fluid ceases to be valid [17]. Therefore, as time progresses, 
vorticity is produced on ever larger scales. It is difficult to study analytically 
the production of vorticity and many works have used dedicated IV-body simu¬ 
lations in order to obtain its magnitude and time/scale dependence [18, 19, 20]. 
In particular, the results of [19, 20] show how vorticity is correlated with the 
density contrast and that it becomes important (as compared to the velocity 
divergence) at small scales and/or at early times. 

The coupling of (12) could then arise if the scalar field is coupled to the vector 
perturbations associated with vorticity. Phenomenologically, the coupling Q will 
be proportional to an invariant contraction of the gradient tensor of the velocity 
field dv.i/dxj (see [18] for a discussion of possible invariants), which, based on 
the results from numerical simulations [19, 20], we can parametrize through a 
power of 8 rn , as in the expression Q = v 8^ above. Note, however, that we are 
parametrizing the nonlinear growth of perturbations (the one correlated with 
the vorticity) with the linear growth function of equation (13). Although this 
is similar in spirit to what is usually done within the spherical top-hat model 
for nonlinear collapse, it is also clear that this approach will only be able to 
describe qualitatively the growth of vorticity. 

The growth function 8 m of equation (13) describes the normalized growth of 
perturbations. Clearly, on smaller scales perturbations are stronger (i.e. the ini¬ 
tial value of 8 m is larger) and nonlinearity is reached at earlier times, as already 
mentioned in the discussion above. It is then natural to question whether a scale 
independent quantity such as S m is able to parametrize the growth of vorticity. 
From the results of [19, Fig. 3] one sees that the magnitude w(r, z) of the vortic¬ 
ity grows at approximately the same rate at different scales. Or, equivalently, at 
the various redshifts examined the scale-dependence of co(r, z) follows a power 
law with approximately constant slope. Therefore, if one were to average the 
vorticity contribution at fixed times, one would obtain an approximately con¬ 
stant factor times the vorticity at some pivot scale r p : (w(r, z)) ss const oj(r p , z). 


3 Not to be confused with the halo spin, see e.g. [16]. 
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The constant factor can then be absorbed into the coupling constant v and, as 
discussed above, uj(r p ,z) can be parametrized using 5 m so that one obtains 
Q = v This should justify the approximation of using the scale-independent 
quantity S m in order to parametrize the scale-dependent onset of vorticity. 

Although we have justified - at least in a first approximation - the use of 
the linear growth function S m as a proxy for the vorticity of the velocity field, it 
would be desirable to quantify directly and more precisely the growth of vorticity 
using numerical simulations. 


4. Dynamical system approach 

As customary [21, 12], let us switch to the following variables: 

V 


x 2 = 


- 


&M 2 X H 2 ’ 


3MpjiJ 2 ’ 


y 2 = 


3 M^H 2 ’ 




3M^H 2 ' 


The Friedmann equation (2) becomes then: 

1 = x 2 + y 2 + + f2 r , 


(14) 


(15) 


so that one variable is not independent. The dark energy density parameter 
and equation of state, and the effective equation of state are, respectively: 


2 2 

_ 2 2 X ~ y 2 2 

= x + y , w^= 2 2 , Weff = X -y + 


3 


(16) 


as: 


x z + y 

We then consider the case of an exponential self-interaction potential such 

V{4>) = V 0 e ~ x<l> / Mp> , (17) 


where A is the dimensionless self-interaction coupling. Using the e-folding time 
variable N = ln(a/ao) = — ln(l + z) (we will denote the derivative with respect 
to N with a prime) equations (3-13) become: 


H' 


1. 


~TF Ax(l + 2 ; — y + ) — 0, 


H 2 
H' 
H 


x 


+ 3 ) x- y 2 = U m , 


^ + ^§Xx)y = 0„ 


H' 


+ ( 2 — + 3 ] — +\/6 v S^ x , 


n' 

r 

X" 

0™ 


2 -gT + 4 ) fi r — 0 , 
H' 


H 


+ 2 ) S' — -Q m 8m — 0 . 


(18) 

(19) 

( 20 ) 
( 21 ) 
( 22 ) 
(23) 
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5. Limiting cases 

We will now discuss the known limiting cases of this model. First, for A = 
v = 0 one recovers the ACDM model, with Hao = 2/(0) = Vo/p c o- 

For v = 0 the system of equations (18-23) reduces to the standard case of 
an uncoupled quintessence field with exponential potential. In particular, if 
A 2 < 3 there is a stable (accelerated for A 2 < 2) attractor with = 1 and 
w<f, = w e ff = — 1 + (scalar field-dominated solution); and, if A 2 > 3, a stable 
attractor with = A and w</, = w e s = 0 (scaling solution). 

For n = 0 one recovers the case of coupled dark energy with Q = v [ 11 , 12 ]. 
In this case, if A 2 < 2 there is an accelerated fixed point with = 1 and 
w,j, = tc e ff = — 1 + This point is stable for A {Q + A) <3. If instead it is 

A(Q + A) > 3, one has a stable fixed point with f w$ = q(q +\)+3 

and w e fi = This fixed point is accelerated if Q < — A or Q > A/2. This 

scaling solution can in principle explain the acceleration of the universe without 
fine tuning. However, it does not produce a viable cosmology as the matter 
dominated era does not set in and matter inhomogeneities do not grow [9]. 

For n yf 0 - the case considered here - it turns out that exists only one fixed 
point, characterized by <$JJ,(oo)fi m (oo) = X/u, for which it is possible to obtain 
also the asymptotic behavior of all the dynamical functions [22] . However, as we 
will see, for a qualitative understanding it is very useful to define “instantaneous” 
fixed points using the instantaneous coupling Q(N) = v5^ n {N). 

6. Evolution 

The basic idea of the model of equations (18-23) is that a coupling to matter 
perturbations makes the scalar field enter a scaling accelerated solution. From 
the discussion above on the limiting cases one concludes that this is achieved 
if the scalar field reaches today the attractor with > w h ere Q = 

jz <5)^(0) = v. At early times, this coupling will be suppressed by the small 
early values of the growth function S m . The growth function provides therefore 
an effectively time-dependent coupling with a step-like functional form. This 
is exactly what is needed in order to have a matter dominated era followed 
by the scaling attractor above, as first advocated by Amendola and Tocchini- 
Valentini [23], in which a transition to the accelerated attractor was obtained 
by means of a (/(-dependent coupling (see also [24]). We will adopt as values for 
A and v the present-day values used in [23]. Within our notation these values 
are {A, v} = {6.5,14.5} and give a scaling accelerated attractor with fl m0 ~ 0.3 
and fl^o — 0.7. Note that, as both couplings are positive, energy is transferred 
from dark energy to dark matter. The model parameters are now fixed (we will 
set n = 2) and only the initial conditions are left to be chosen. 

Figure 1 (left column) shows the dynamics for the initial conditions Xi = 0, 
Di = 10 -5 - which correspond to = 10~ 1Q - and H r ,; = 0.998 at the initial 
redshift log 10 (l + z) = 6. Figure 1 (right column) shows instead the dynamics 
for Xi = 0, yi = 10 -1 ’ 5 - which correspond to = 10 -3 - and Q ri = 0.998 
at the same initial redshift. The initial conditions for equation (23) have to be 
given taking into account the early radiation dominated era. Specifically, they 
are 

~Sm(N i )= e ^ + l e ^, 5' m (N i ) = e N *, (24) 
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Figure 1: Comparison of the dynamical evolution of the model presented in this paper with 
parameters {A, v, n} = {6.5,14.5, 2} for two very different sets of initial conditions. The plots 
of the left column are relative to initial conditions {f^,Q r i} = {10 10 ,0.998} at the initial 
redshift log 10 (l + z) =6, while the plots of the right column are relative to {fi,£j,Q r j} — 
{10 -3 ,0.998} at the same initial redshift. First row: trajectory in the {x,y} phase space of 
the solution of equations (18-23) (solid blue) and of the instantaneous attractor (dotted red, 
see section 5). The black curve is the border of the allowed region given by y = \/l — x 2 . 
Note how the model converges to the present-day accelerated scaling instantaneous attractor. 
Second row: redshift evolution of the relevant density parameters. The dynamics features the 
expected series of radiation, matter and dark-energy dominated eras. Third row: evolution 
of the dark energy equation of state and of the effective equation of state. Also shown for 
comparison is the effective equation of state for a flat ACDM model with the same present- 
day values of the density parameters. Fourth row: evolution of the growth rate f g = S'^/Sm 
compared to the growth rate relative to the corresponding flat ACDM model defined above. 
As can be seen by comparing the plots of the left column with those of the right column the 
model of this paper does not suffer from the fine tuning of initial conditions: the dynamics is 
basically unaltered if the initial dark energy density is changed by seven orders of magnitude. 
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where N eq is the e-folding time at which Q r (N eq ) = Cl m (N eq ) [for details, see 
25]. 

In figure 1 (first row, left panel), the model trajectory (solid blue) in the 
{x,y} phase space reaches the point {xo,yo} = {0.0468,0.833}, the instanta¬ 
neous attractor trajectory (dotted red) reaches instead the point {0.0583, 0.833}. 
The field is basically on the instantaneous attractor, which is then useful for un¬ 
derstanding the dynamics. Also, it is clear from the plot how the initial attractor 
was frozen (and not accelerated) for a long time allowing the model trajectory 
to spiral around it. Then the effective coupling started making the attractor 
move towards its present location. 

If the initial condition for is changed to the substantially different value 
= Ui = 10 -3 = f l mi - at this very step one usually encounters the 
fine-tuning problem as shown in figure 2 - one has basically the same tra¬ 
jectory, as shown in figure 1 (first row, right panel). The end point is now at 
{0.0467, 0.834} ~ {xq, yo}, showing how insensitive the evolution is with respect 
to the initial conditions. The same conclusion holds for any viable values of A 
and v [22]. 

The second row of figure 1 shows the redshift evolution of the density pa¬ 
rameters, and that this model achieves the matter dominated epoch after the 
radiation dominated one. After the matter dominated era, the dark energy 
dominated era sets in as the scaling attractor with 0 ~ 0.3 and — 0.7 
is reached. The third row of figure 1 shows the redshift evolution of the dark 
energy equation of state and of the effective equation of state. Also shown is 
the effective equation of state for a flat ACDM model with the same present- 
day values of the density parameters. The comparison of the plots of the left 
column with those of the right column shows that this model does not suffer 
from the fine tuning of initial conditions: the dynamics is basically unaltered 
- see for example the evolution of the effective equation of state - if the initial 
dark energy density is changed by seven orders of magnitude. 

The dynamics shown in the second and third rows of figure 1 is qualitatively 
similar to the standard coupled dark energy evolution if the parameters are 
chosen in order to approach an accelerated fixed point with flj, = 1 and w $ = 
w e g = — 1 + Aj-. However, the latter model suffers severely from the coincidence 
problem and its related fine tuning of initial conditions. For the ACDM model 
and standard coupled dark energy models to work, one has to set the initial 
density parameter precisely to a value of about ~ 10 _2 °, tens of orders of 
magnitude smaller than the initial energy density of matter and radiation. This 
is shown in figure 2 where the evolution of the density parameters for {A, u, n} = 
{0.1, 0.3,0} is shown for different initial values of the dark energy density. 

Finally, the bottom row of figure 1 shows the redshift evolution of the growth 
rate f g = S^/Sm compared to the growth rate / g ,a relative to the corresponding 
flat ACDM model defined above. As it can be seen, the overall redshift depen¬ 
dence of the growth rate is qualitatively similar to the ACDM growth rate, even 
though moderately suppressed due to the fact that - for this specific choice of 
parameters - dark energy was subdominant but not negligible in the past. Non 
negligible dark energy in the past could affect the CMB as Early Dark Energy 
studies have shown [26]. 







fi 01 =io - 22 V 

n r 


. X 2 

-- --- 
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0 1 “““ 2 3 4 5 6 


logio(z+l) 


Figure 2: Redshift evolution of the density parameters for the standard coupled dark energy 
model (n = 0 within the notation of this paper) with {A, v} = {0.1, 0.3}. The accelerated 
fixed point has = 1 and = w e R ~ —0.997. The initial conditions are = 10 -20 
(top), City = 10 -21 (middle), = 10 —22 (bottm) and = 0.998 at the initial redshift 
logio(l + 2 ;) = 6. It is clear that the middle plot features the desired series of radiation, matter 
and dark-energy dominated epochs, though at the price of a severe fine tuning of the initial 
conditions. On the contrary, a similar evolution can be obtained by the model presented in 
this paper without tuning the initial conditions, see figure 1. 
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7. Conclusions 


We have discussed the possibility that dark energy could effectively couple 
to matter inhomogeneities. Although only phenomenological, this model is in¬ 
teresting as it is possible to obtain a viable cosmology free of the coincidence 
problem for a wide range of parameters. In particular, the current density of 
dark energy is given by the value of the coupling parameters A and v and is 
insensitive to the precise value of the initial conditions for the scalar field. 

We proposed that this phenomenological interaction is due to the coupling 
of the scalar field to the vector perturbations associated with vorticity. Vorticity 
is indeed produced at very non-linear stages of the dark matter collapse when 
the single-stream irrotational description of the dark matter fluid ceases to be 
valid. Therefore, vorticity is highly correlated with matter perturbations, as 
also seen in dedicated numerical simulations which could be used to precisely 
evaluate the growth of the average vorticity in the universe. Also left to future 
investigation is a comparison with present and future cosmological data in order 
to see if the model is viable and, in case, if upcoming observations could rule it 
out. 
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